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Abstract 

The observed dark matter phenomenon is attributed to the presence 
of a gas of wormholes. We show that due to topological polarization ef- 
fects the background density of baryons generates non- vanishing values for 
wormhole rest masses. We infer basic formulas for the scattering section 
between baryons and wormholes and equations of motion. Such equations 
are then used for the kinetic and hydrodynamic description of the gas of 
wormholes. In the Newtonian approximation we consider the behavior of 
density perturbations and show that at very large distances wormholes 
behave exactly like heavy non-baryon particles, thus reproducing all fea- 
tures of CDM models. At smaller scales (at galaxies) wormholes strongly 
interact with baryons and cure the problem of cusps. We also show that 
collisions of wormholes and baryons lead to some additional damping of 
the Jeans instability in baryons. 



1 Introduction 

The nature of Dark Matter (DM) represents one of the most important and yet 
unsolved problems of the modern astrophysics. Indeed, while the presence of 
DM has long been known [I] and represents a well established fact (e.g., see 
[H 15 and references therein) , there is no common agreement about what DM is. 
In the simplest picture DM represents some non-baryonic particles (predicted 
numerously by particle physics) which should be sufficiently heavy to be cold at 
the moment of recombination and those give the basis to the standard (cold dark 
matter) CDM models. The latter turn out to be very successful in reproducing 
properties of the Universe at very large scales (where perturbations are still 
on the linear stage of the development) which led to a wide-spread optimistic 
believe that non-baryonic particles provide indeed an adequate content of DM. 

However the success of CDM models at very large scales is accompanied with 
a failure at smaller (of the galaxies size) scales. Indeed, cold particles which 
interact only by gravity should necessary form cusps {poM ^ ^/r) in centers 
of galaxieq^ [4] (see also [5] where the problem of cusps in CDM is discussed 



^The presence of cusps formed by the development of adiabatic perturbations follows 
straightforwardly from the conservation of the circulation theorem in the hydrodynamics. By 
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in more detail), while observations [B] definitely show the cored [poM ^ const) 
distribution. The only way to destroy the cusp and get the cored distribution 
is to introduce some self-interaction in DM or to consider warm DM. Both 
possibilities are rejected at large scales by observing AT/T spectrum (e.g., see 
[3] and references therein) . By other words DM displays so non-trivial properties 
(it is warm or self-interacting in galaxies, however it was cold at the moment 
of recombination and it is still cold on larger (than galaxies) scales) that it is 
difficult to find particles capable of reconciling such observations. 

These facts support the constant interest to different alternatives of the 
DM hypothesis which interpret the observed discrepancy between luminous and 
gravitational masses as a violation of the law of gravity. Such violations (or 
modifications of general relativity (GR)) have widely been discussed, e.g., see 
[ll [8]. However, it turns out to be rather difficult to get a modification of GR 
which is flexible enough to reconcile all the variety of the observed DM halos. 
Moreover, the weak lensing observations of a cluster merge [HI seem to reject 
most of modifications of GR in which a non-standard gravity force scales with 
baryonic mass. 

The more viable picture of DM phenomena was suggested by [10] (see also 
references therein) and developed recently by [13] [TB]. It is based on the fact 
that on the very early (quantum) stage the Universe should have a foam-like 
topological structure [T7] . There are no convincing theoretical arguments of why 
such a foamed structure should decay upon the quantum stage - relics of the 
quantum stage foam might very well survive the cosmological expansion, thus 
creating a certain distribution of wormholes in the Friedman space. Moreover, 
the infiationary stage in the past [18] should enormously stretch characteristic 
scales of the relic foam. The foam-like structure, in turn, was shown to be flex- 
ible enough to account for the all the variety of DM phenomena [THI [Ilj ; for 
parameters of the foam may arbitrary vary in space to produce the observed 
variety of DM halos in galaxies (e.g., the universal rotation curve for spirals con- 
structed by [13] for the foamed Universe perfectly fits observations). Moreover, 
the topological origin of DM phenomena means that the DM halos surrounding 
point-like sources appear due to the scattering on topological defects and if a 
source radiates, such a halo turns out to be luminous too [TB] which seems to 
be the only way to explain naturally the observed absence of DM fraction in 
intracluster gas clouds |9j. 

The foam-like structure of the Universe is represented by the gas of worm- 
holes randomly distributed in space. It was demonstrated recently by [Tl] that 
in the presence of such a gas every point source turns out to be surrounded with a 
" dark halo" which possesses both signs depending on scales and the background 
distribution of wormholes. Therefore, it seems to be not quite clear for readers 
whether wormholes produce the necessary CDM picture, or they merely suggest 
some additional effects. Moreover, there still exists some wide-spread mistaken 
opinion that wormholes lead to non-Gaussian perturbations in the analogy with 

other words the fact that the distribution of DM should have cusps in galaxies is equivalent 
to the fact that DM should represent cold non-baryonic particles. 
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topological defects of another kind (strings, monopoles, etc.). In the present pa- 
per we clarify such problems by means of considering the development of density 
perturbations in the gas of wormholes. We demonstrate that at very large scales 
wormholes behave exactly like very heavy particles and thus reproducing all the 
predictions of CDM. They however predict an additional specific damping in 
the development of baryon perturbations. Moreover, at smaller scales the non- 
linear stage of the evolution of perturbations essentially diverges from that in 
CDM. At small scales there exists a rather strong non-gravitational interaction 
between particles and wormholes (due to the mutual scattering) which surely 
cure the problem of cusps in galaxies. 

The complete analysis of the gravitational dynamics of wormholes is rather 
complicated. Therefore, in the present paper we restrict ourself with the Newton 
approximation to derive basic equations which govern the dynamics of baryons 
and wormholes. While the problem of the generalization to the General Rela- 
tivity we leave for the future research. 

The paper is organized as follows. In Sec. 2 we introduce a wormhole and 
describe it's general properties. In Sec. 3 we show that in the Friedman model 
wormholes acquire non-vanishing rest masses. In Sec. 4 we describe the process 
of the scattering of particles on a wormhole. In Sec. 5 we infer the motion 
equations for particles and wormholes in the Newtonian approximation for the 
expanding reference system. In Sec. 6 we introduce the system of Boltzmann - 
Vlasov equations which describes kinetics of particles, wormholes, and collisions. 
In Sec. 7 we infer the non-relativistic hydrodynamic equations with corrections 
for the collisions between particles and wormholes. In Sec. 8 we calculate basic 
kinetic coefficients which describe the collisions. In Sec. 9 we consider the be- 
havior of linear perturbations in the Newtonian approximation. In particular, 
we explicitly demonstrate that at large distances wormholes behaves exactly 
like heavy non-baryon particles reproducing thus the standard CDM picture. 
While the strong coupling at smaller scales cures the basic failure of CDM (i.e., 
it removes cusps in galaxies). In the last section we discuss results obtained and 
show further perspectives. 



2 Wormholes 

The simplest wormhole is described by the metric 

ds"^ = (?dt'^ - (r) Sc,^dx"dx^, (1) 

where 

h{r) = l + e{b-r)(^^-l^ (2) 

and 9 (x) is the step function (0 (x) = as .t < and 9 (x) = 1 as .x > 
). We point out the obvious properties h'^ = 9{h~r) and /i"^ = 

5 (6 — r) ( ^ j {b — r)' + 9 [b — r) i^j . Such a wormhole has vanishing throat 
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length. Indeed, for the region r > b h = 1 and the metric coincides merely with 
that in the Minkowsky space, while the region r < b, upon the obvious trans- 
formations — ^x", gives the same region y > b with the same flat metric 
ds^ — c^dt^ — 5ai3dy°'dy^ . Therefore, both regions r > b and r < b represent 
similar portions of two Minkowsky spaces glued at the surface of the sphere 
with the center at the origin r = and the radius r — b. Thus, formally, such a 
space can be described with the ordinary double- valued flat metric in the region 
r± > 6 as 

ds^ = c^dt^ - Scpdxidxi, (3) 

where the sign ± in coordinates x^. stands to describe two different sheets of 
space. 

A generalization appears when we change the step function 6 (x) with any 
smooth function 6 (x) which has the same property 9 (x) — )■ as x ^ — oo and 
9 (x) — >■ 1 as a; — > 0. On the contrary to the previous case such a wormhole will 
have a non- vanishing throat length. However in the last case the consideration 
differs in details, while general features remain the same. We also point out 
that such wormholes have vanishing mass. In general, one may also insert a 
non- vanishing mass to the wormhole as it is described in 19,. 

First of all we consider the stress energy tensor which produces such a worm- 
hole. It can be found from the Einstein equation = — ^6^R. Since the 
metric H]) does not depend on time we flnd 

i?° = i?o=0, R^^ = -P^ = -'isib-r){7^^n^-Si} (4) 

where n" = ric = x°' /r is the outer normal to the sphere and therefore 

T;:^-ls{b-r){6;~A;) (5) 

where A[] = A° = 0, and A^^ = - naU^ = YjA=i,2 (^t^"^^ where are 
two unite tangent to the sphere vectors e^n" — 0. Thus, the effective source 
can be considered as a mixture of a negative " cosmological constant" T^'^ = 
— 15 (6 — r) (5^ and two components of a perfect fluid T^^" = — 1(5 (6 — r) u^u^'^ 
with zero pressure p = 0, energy density e — — 1(5(5 — r), and the velocity 
= (O, e^'^). The velocities have purely space-like (tangent to the surface of 
the sphere) components u^w^^ = — 1. This matter is concentrated only on the 
surface of the sphere due to the multiplier 5 (b — r). All such sources represent 
an exotic matter which cannot be constructed from actual particles. Recall, 
that all real particles have time-like velocities u^u^ = 1 or isotropic u^u^ = 
(if the rest mass vanishes m = 0). However, this property is violated for virtual 
particles. This gives us a hint that vacuum polarization effects can in principle 
be collected to organize such a form of matter. By other words to construct a 
wormhole we have to organize a negative cosmological term on the surface 
and then to rotate the surface with the space-like velocity in the two different 
directions. From the other hand the presence of a wormhole leads to the vacuum 
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polarization and therefore to the exotic matter (e.g., see [lOlISI])- Here we leave 
aside the important problem of what precedes the egg (wormhole) or the hen 
(exotic matter). 

The inner and outer regions of the sphere are equal, which gives the 
possibility to construct a wormhole which connects regions in the same space 
(instead of two independent spaces). This is achieved by the the identification 
(gluing) in (jSj of the two spaces by means of the use of possible motions of the 
flat space. Let R+ be the position of the sphere in coordinates x", then the 
gluing is the rule 

xl^R1 + C/^ (a;^ - i?^) (6) 

where Up G 0(3), which represents the composition of a translation and a 
rotation of the flat space. In terms of common coordinates such a wormhole 
represents the standard flat space in which the two spheres S^. (with centers at 
positions i?^) are glued by the rule ^ 

If we neglect the higher order images with respect to the transformations 
Q = (which are the reflections with respect to spheres 5'^), then we can 

use the metric ([l} in which we merely replace h with the function 

h (r, i?±) = 1 + (6 - l^+l) (^^-?j+e{b- le-l) (^^ - l) (7) 

where = x" — R% and ^'^ ~ (J7^)*^^£. In this case the physically admissible 
region of space is the outer region of the two spheres S'^. , while the inner regions 
represent only additional images of the outer space and have no meaningj. 

Thus, we see that the simplest wormhole may be described by a set of pa- 
rameters r] = i?" , b, . We point out that in the general case a worm- 
hole possesses a rather complex properties, e.g., parameters possess a dynamics 
rj = rj{t) and, moreover, (6, U^) are functions of coordinates x°' which describe 
the specific structure of the wormhole. 

The generalization to a set of wormholes is given by 

/i(r,i?±) = l+ HK~\£,na\){^-l) (8) 

and the identification of points at joint spheres S^^^ as ^"_|_ = C/^o^f_. 



3 The rest mass of a wormhole in the expanding 
Universe 

As it was demonstrated previously (e.g., see for details [11]) the presence of a 
gas of wormholes leads to a specific topological polarizability of space. This can 

■^If we do not neglect the higher order images, then within every sphere will appear a count- 
able set of spheres which represent additional images (copies) of the outer (physically 
admissible) region of space. 
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be described as a bias of point-like sources of gravity 



6{r ~ ro) 6{r ~ ro) +B{r,ro), 

where the bias gives the polarization mass generated on throats. The bias 
consists of two terms B = Bq + Bi, where Bq resembles the bias of spherical 
mirrors and gives the positive contribution (anti-screening) to the total mass of 
a particle 

Bo{r)^ J2 ^\s{r-rn,.)-S{r-Rr,,,)\, (9) 

n,(T—± 

where r„ + — ± + ,_ °^ ,o U^^(fo — Rn while the rest part gives pure 
screening and is given by 

B, (r) = ^" {j^ J^) ['^(^^- ^^".-) - '^(^^- ^^".+)] ■ (10) 

In the homogeneous Universe (i.e., in the case of a homogeneous distribution of 
matter) the second part disappears automatically, while the first part generates 
the rest mass for every wormhole. Indeed, to demonstrate this we rewrite the 
bias ([9]) in the form (see for details [TT| ) 



where 



'R^ 



<FR, 



Mr) = 1 ^ [Kp (r, r) + [r, f 

(R+R-) - J b^U^^F (i?±, 6, U) dadU, 

and F (ry) is the distribution of wormholes in the configuration space. In the 
case of the homogeneous and isotropic distribution of wormholes, averaging over 
the rotation matrix gives H^^ (i?+^_R_) = ^Sai3(f>{\R+ ~ R-\), where (l>{X) = 
J b^F {X, b) db. Therefore, the bias reduces to the form 

47r — 

Boir,ro)^—nb^{r~fo) (12) 

where n is the density of throats, while ^nb^ is the portion of the unit volume 
which is cut by throats. In the above formula the only difference from spherical 
mirrors (the multiplier 1/3) appears due to averaging over the rotation matrix. 
By other words every wormhole increases the mass of a point-like source in 
proportion to the volume which the wormhole cuts from the space. When the 
space is filled with homogeneously distributed matter, every wormhole acquires 
the rest mass Al^ [b] = ^b^ p [p is the mean density). 
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All the above consideration can be omitted if we note that the wormhole 
possesses straightforward generalization to the case of the expanding Universe. 
Indeed, the metric is 

ds2 ^ „2 ^^^2 _ ^2 S^^dx^dx/^) . (13) 

In this case the wormhole expands with the space and acquires a non- vanishing 
rest mass 

M„ = -ira^b^pb (14) 

where pi, is the mean density of matter in the Universe. By other words the 
gluing procedure cuts two equal homogeneous portions of space (the two spheres 
Sj.) filled with a homogeneously distributed matter. Thus, the total amount 
of matter diminishes, while spheres acquire the rest masses which compensate 
the distribution of the effective density of matter to the homogeneitjo. We 
point out that the metric (fT3|) . with ([8]) taken into account, directly shows 
that an arbitrary primordial distribution of worniholes in space agrees with 
the visible homogeneity of space (the Friedman model) as it was first stated 
by [To]. An inhomogeneous distribution of wormholes in space (e.g., fractal 
distribution) leads to the dual fractal distribution of particles (for particles may 
occupy only the physically admissible region of space), while the total effective 
density remains perfectly homogeneous. 

From the astrophysical standpoint such a mass is much smaller than that 
of a typical object (e.g., if the throat radius is rth = ab ^ i?o, and pb = Pcr, 
where i?o is the Sun radius and pcr is the critical density, then the mass has 
the order M^, ^ lOkg). Thus, gravitational effects of such objects in Solar sys- 
tems is merely negligible. However in the dynamics of larger systems (galaxies, 
clusters, etc.) they become the more and more visible , e.g., as the presence of 
dark matter. We point out that such dark particles (wormholes) are extremely 
heavy from the particle physics point of view. It is clear that the picture where 
the Universe is filled with a gas of " particles" whose rest masses have the order 
~ 10kg perfectly fits the basic Cold Dark Matter model (CDM) widely accepted 
today. And as it is well known, predictions of CDM perfectly agree with ob- 
servations at very large (clusters of galaxies and higher) scales. However, while 
the standard CDM (heavy non-baryon particles) makes a wrong prediction at 
smaller scales (e.g., it predicts cusps pdm ^ ^/r in centers of galaxies, when 
observations demonstrate the cored pdm const distribution), we may expect 
that wormholes will cure such a problem. Indeed, since the two conjugated 
spheres Sj_ (throats) represent merely the " same" region of space we may state 
that the local density at S'^ coincides exactly with that at S"^. Let do be the 
typical distance between throats. Thus, if one throat gets into the central region 
of a galaxy while the rest throat is sufficiently far from the center, then the to- 
tal density will be somewhat smoothed. And the typical scale of the smoothing 

^We also point out that rotating wormholes will bo described by the tensor of inertia 
= IS^is, with / = §a2fe2M„. 
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(i.e. the minimal core radius) will be of the order of do- Of course, the rigor- 
ous consideration of this problem requires considering the proper dynamics of 
wormholes which we start in the next section. 



4 The collision with a wormhole 

Consider the simplest construction of a wormhole in space as follows. Let us 
fold the space at the plane z — 0. Then we can use the simplest metric 
which connects the two half-spaces z > and z < 0. Let the velocity of the 
wormhole is V which, due to the above symmetry (i.e., Z+ — ~Z^, where 
B.± = {X±, Y±, Z±) are the positions of throats), corresponds to the case when 
one throat (say S'^) moves with the velocity V- = V , while the other S\ moves 
with V+ = {Vx,Vy, — Vz). For the sake of simplicity we consider everywhere the 
non-relativistic case (i.e., V c). In the case of a general wormhole the relation 
between V+ and V- is given by the same relation ([6|) 

= Uj^vt (15) 

By other words since both spheres represent the same region of space their 
velocities are rigidly connected (the visible doubling of the number of degrees 
of freedom related to wormholes is fictitious). 

Consider an incident on particle with the rest mass m and the initial 
velocity v. Then the scattering of the particle leads to the transformatiorQ 

n° ^ = R± ^ R'^ ^ R± (16) 

where = (x" — i?^) /b are points at spheres Sj., and the transformation of 
velocities 

Here (un) — (ufi) denotes the ordinary scalar product and 



and the rotational matrix Up € 0(3) corresponds to the above symmetry (i.e., 
the reflection with respect to the plane z = Q, z' — —z). 

*This law of the transformation can be easily obtained as follows. We notice that the 
metric {T} merely connects two equal flat spaces. Then the scattering of a particle merely 
coinsides with the elastic scattering on a solid ball with simultaneous transport of the particle 
from one sheet of space to the other sheet. 
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We point out that the above law of transformation ([T5 1) - (|18p holds in the 
most general case (the case of an arbitrary gluing, i.e., an arbitrary Up S 0(3)) 
as well. In the limit — oo this law has been considered by us in |15) . 

It is easy to see that the law of transformation P^ - ((T5| conserves the total 
energy 

mv^ + MV^ = mv'^ + MV_^ 

where the kinetic energy of the wormhole is accounted only for once. It maybe 
convenient to consider throats as an independent objects. Then, due to the rigid 
relation we have to introduce the bound energy e = —^MV^ = —^MVf 
which gives the conservation of the total energy in the form 



mv^ MVl MVl _ MV^ MV^ 



(19) 



It is important that if we forget about the factorization (i.e., the identi- 
fication of points or the gluing) (|6]), then the above transformation does not 
conserve the momentum. Instead we get the conservation in the form 

+ mv'" = {MV^ + mv°') . (20) 

It happens due to the fact that the space with the wormhole does not possess 
the translation symmetry. It possesses however the somewhat entangled (with 
the gluing ([6]) taken into account) symmetry. 

We point out that when we consider dynamics of wormholes the matrix 
becomes function on time (formation of gravitationally bounded objects leads 
naturally to origin of rotational motions). 



5 Equations of motion 

As it was shown above a general wormhole possesses a rest mass and, therefore, 
it should move in space as an ordinary test particle. For the sake of simplicity 
we shall use an approximation when wormholes can be considered as point-like 
objects, i.e., we neglect the size of wormholes. For cosmological applications this 
represents a rather good approximation, while for some astrophysical problems 
it can be not sufficient. 

The presence of a gas wormholes in the Universe (i.e. of the complexity of the 
topological structure of space) leads to an enormous complexity in consideration 
of the dynamics (due to boundary conditions at wormholes). To simplify the 
problem we shall use the Newton's equations, while the generalization to the 
relativistic case we leave for the future. Moreover, as it is well known (e.g., see 
[22] ) the Newton's approach is acceptable for a rather huge range of scales. It 
will be also convenient to introduce the expanding reference system from the 
very beginning. 
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5.1 Newton's equation 

When the gravitational field is rather weak, the Einstein equations reduce to 

Wl<S> ^ AttG ( p + ^] (21) 



where « Rqq is the time component of the Ricci tensor. The equations of 
motion for a test particle are 

which have the range of the applicability w << c, or (equivalently) at distances 
R << cH~^{H = d/a is the Hubble constant). The lower boundary is given by 
the gravitational radius for compact objects (black holes). 

Consider the homogeneous component P{t) and p(t) and the homogeneous 
expansion of space f= a{t)x. Then the equation (|2ip gives 



c? ) 



while the equations of motion ([22]) lead to the standard cosmological equation 

Here the point a; = corresponds to the observer (which gives ~ r^) and the 
applicability is restricted by 

5oo = l + 2$/c2 2$<<c2 

which means that distances r cannot be too large. 

5.2 Peculiar velocity 

Consider now the expanding reference system r — a(t)x. Then the total velocity 
is 

u = ax + xa = ax + Hr 

where ax is the peculiar velocity, while Hf gives the standard Hubble expansion. 
For a particle the equations of motion can be obtained from the Lagrangian 

m) 

C = —m yax + xdj — ni^{x, t). 

Let us use the canonical transformation 

r . r ^'^ I 1-2 
L ^ L — , w = -maax 

dt 2 
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which transforms the Lagrangian function to the form 



where 



r ^ 2 '-2 



1 2 
(/3 = $ + -aaa; 



and new Newtonian potential <y9 satisfy the equation 



\ y a 

where the gradient is taken over x. Then by means of use of the cosmological 
equation ([23l) it transforms the Poisson equation to the form 

where (5p = p(x^ t) — Pb(t). The equations of motions for a test particle become 

2U dp 

p = ma X, ~ —nriS/Lp (24) 

where the peculiar velocity v = ax relates to the momentum as v = p/ma and 
by virtue of (p4|) obeys the equation 

dv a ^ 1 _± 
— - + = — V(p. 
at a a 

E.g., if V(/3 = 0, then p = const and the velocity changes as v ^ — , i.e., 
^l/ait). 

Thus, to describe the cosmological model filled with a set of particles and 
wormholes we get the system of equations as follows: 

4 



■d = --TiG ypb{t) + 

= 4ttG ( dp{x, t) + ) (25) 



mAa^XA = PA 
PA = -mAV^p{xA) 

where A = 1, ...,N numerates particles or wormholes. 

This system of equations looks as the ordinary one but we should take into 
account that wormholes insert some complex (depending on time) identification 
of points (boundary conditions for and xa)- Moreover, since wormholes are 
described by a larger number of parameters than a particle, i.e., 77^, "HA {11 = 
(^R+i R- , b, Vjf^ ) to complete this system we should add equations which govern 

the evolution of the rest parameters {bA and Up a)- For the sake of simplicity 
in the present paper we shall neglect the possible evolution of h and Up . 
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6 Boltzmann -Vlasov type equations 

In what follows we shall use the standard methods (e.g., see [23]). Let us 
introduce the density of particles in the phase space as r„j = (x, p) , 

and the number density of wormholes 

dN.u, (7) = Eu,{T^,j,t)dTw. 

here, since we assume the part of parameters 7 = (6, Up) be fixed (non-dynamical 

ones), dTw = drjdpri, — (rj^prj), and, 77 = (i?-|_,i?_). We also note that ac- 
cording to (jl5p only one momentum is free. It is convenient to consider both 
momenta Prj = (P+,P-) as independent, while the relation is accounted 
for by a delta- function - type multiplier in . Then the matter density can be 
expressed via fm{x,p,t) as 

p^{x,t) = — JS{x~x')f„^{T'^,t)dr'„, = pb{t)[l + S^{x,t)] (26) 

and analogously the contribution from wormholes is (here M (7) — ^Trb^a^pi,) 

Pw{x,t) = X! y {x-R„^ F.u]{T^,'^,t)dr^d'y. 

We point out that in the proper (expanding) reference system the mean number 
of particles in the volume element (Ax) remains constant, i.e., a^pb ^ const 
and so do the masses of wormholes M (7) = M (b) ^ const. 

In Fw{Tw, 7, t) indexes w and 7 numerate merely the sort of wormholes. It is 
convenient to introduce the reduced distribution function for wormhole throats 
as follows 

U(x,P,-f,t] = Y. 1 5 (x - Ra) 5 (P - Pj\ F^{T^,t)dT^ 

which defines the mass density in the form analogous that of particles 
M(7) 



p^{x,t)^ -nlf^^(^S,P,-f,t')d^Pd^^Pb{t)[l + 5^{x,t)]. (27) 



Then the system of Boltzmann- Vlasov equations are eqs. (|26| . p7| 

dfA ^ P 
dt niAOi 



1^ + ;;r:7;2 v/a ~ '^av^— = ^ stfAs, (28) 



^ B 

2 

and the generalized Poisson equation (we set SP/c^ = -^Sp ^ Sp otherwise one 
has to add the coefficient (1 -I- Sc^/c^) before Sp) 

^A(p = AttG (^Sp{x, t)+ J B {x, x') Sp'd^x'^ . (29) 
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In the last equation the bias B {x, x') accounts for the topological polarizability 
of space in the presence of the gas of wormholes (i.e., the proper boundary 
conditions at wormhole throats). As it was shown by [TT] it expresses completely 
via Fw{Tyj,^,t) and is given by ©, (ITU|) . In (^5)) the index A denotes either m 
in the case of particles, or (if, 7) in the case of wormhole throats. 

Rigorously speaking the above system is not complete, since the collisions 
between particles and a wormhole involve both throats 

StfAB ^ Jw {T'a, T'b.Ta, Tb) [J'aF'b - IaFb] dV^dT'sdYB (30) 

where the scattering matrix w (F^, F^, Fa, T^) stands for the transformation 
law (Fa,Fb) ^ (r^^r's) which is defined by ([HI)- (El)- The function Fa obeys 
the equation similar to the two-point correlation function 

^ + E {^^'^ " ^aV.^^) Fa^J2 '^P^b. (31) 

Then equations (^5)1 are obtained by the integration of the last equation over 
the extra variables. 



7 Fluid equations 

Upon integrating (|28|) over / d^p we get 



dpAa^ , l.vf^,^,).^^, (32) 



dt rriA a \ rriA 

which represents the first hydrodynamic equation (the discontinuity equation)). 
The term Da ~ J stfAd^p we define latter on, while we used the notions for 
the density 

fAd^p^P^^^^ 

ruA 

the peculiar macroscopic velocity u{x, t) 

- f--^\ 1 JpfAd^p 

UA(x,t) = (33) 

ruAa J jAd-^p 

and pA{x,t) — pi,{t){l + SA{x,t)). Due to the property pi,a^ = const this 
equation can be rewritten as 



^SAi£,t) + -Vp(il + SA)u^A)=^DA. (34) 
ot a \ / pbtt'^ 



We point out that in the standard hydrodynamics the terms of the type Da 
vanish due to the local conservation laws (of the number of particles, momentum, 
energy). In collisions of particles and wormholes such terms retain. During the 
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scattering of particles on wormholes the number of throats conserves locally 
which gives D^, = 0. Though, in general, there remain processes of throats 
merging which produce some value < 0. It is also easy to see that all such 
terms vanish in the absence of perturbations (for the homogeneous background). 

Multiplying the equation (|28p with pa and integrating over momenta d^p we 
get the second (Eiller) equation 

{a^PAul) + pAa^'Vo.V + a'V^ f<w>A) + a'V„P = (35) 



dt 

where 

Qa = / p"stfAd'p 
and we used the obvious definitions 



.0 ^lIa „/a„2\ _„"„/3 , Pa 



3 PA 



(36) 



where Pa is the pressure. The system of equations (|32|) and (|35|) represents 
the basic hydrodynamic equations. In the same way one can also add to the 
above system the equation for the heat transport (i.e., the next momentum for 
the internal energy to^(5w^)/2). For the sake of simplicity we do not consider 
the heat equation here. For adiabatic linear perturbations such equation is not 
important. It however becomes important more latter on non-linear stages of 
the formation of astrophysical objects. 

In view of phc? = const ([35]) maybe rewritten as {H = d/a is the Hubble 
constant) 

^+h] ((1 + 5a) ul) + -{! + 5a) V„¥'+ 



dt y " a 

+ -Wfi ((1 + 5a) + — VoPa = -^Ql- 

a \ / apt a Pb 

Let us take the divergency from the last equation and in view of (j34p we find 
the master equations which govern the evolution of density perturbations (here 
Qa = Vo.Q'X) 

d ^ ^\ f d ^ m „ \ 1 
dt J \dt pba 



-.+^h){-5a- —Da = -2 V„ (1 + 5a) V^^^ 



+— Va ((1 + 5a) + -^^^Pa ~ -^Qa. 

Retaining just linear terms we get the equation for linear perturbations as 
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where we used the relation 5P = c^Sp and Cg is the sound speed. We also have 
to add the Poisson equation 




Thus, the system ((37|) and (l38l) represents the master equation which de- 
scribes the development of linear adiabatic perturbations in the presence of the 
gas of wormholes. This system represents the standard system for CDM model 
where wormholes play the role of dark matter particles. The difference appears 
however due to the additional terms Da, Qa, and the bias B(x,x'). For ho- 
mogeneous background some terms disappear. As we shall see latter on they 
disappear in the long- wave limit as well i.e., as fc ^ fcj (where kj is the Jeans 
wavelength whose definition in the presence of wormholes somewhat differs from 
the standard one). 



8 Kinetic coefficients D and Q 

Consider an arbitrary point x on the sphere i.e., x e S- and therefore = 

^x — R-^ ~ The gluing procedure transforms this point into a conjugated 

point x' £ 5*+ which has the form x' = + ^+ where ^+ relates to ^_ by some 
rotation ^" = Vj^Sj^. Let r„ = (i?±, bw, U, V±) denote the set of parameters of 
the wormhole. The scattering matrix W is presented a,s W = W+ + W- where 
W± corresponds to which of throats absorbs particles. Due to the obvious 
symmetry between throats this gives merely the factor 2 in final expressions. 
Then we find from W±iTm,r^,T'^,T'^,) = |u±| (7±(r,„, T^, T^^, T^, ) 

where ^ 

a_ = si^+- b^) 6 (r;, - Tra,-)S (r'„ - r^) - (39) 

-S{^^-b^)S{T:,,-T„,)S{r'^^T.^) 

and analogous term (t+ with the obvious replacement ( > + and f/ — > U~^), 

where we used the notions as follows ^± = ^ — R±, n± ~ C±/&ii>, 

r„j,_ = (a;_,p_) , = {R±,Vl,b^,U) , 

and the relations 

x"_^R'^ + U^f^, 
p_^U[p-— — ■ (u_n_)n_ , (40) 

Vi^V±+ {u±n±) n±, (41) 

Mw + m 

and 
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For astrophysical needs wormholes may be considered as point-like objects. 
This is achieved by the replacement in ([39)) 



Si^±-b^)^ AnblS - xj (42) 

where S (x) is the 3-dimensional delta function J S {x) (Px — 1. Thus the above 
expressions completely define the scattering matrix W and the scattering term 

stfAB ^ Jw {T'a, Vb^Ta^Tb) [I'aF'b - IaFb] dT'j,dV'BdVB. (43) 
8.1 The scattering of particles 

Consider first the terms Dm and Qm- Since we suppose that A/^, ^ m, it fol- 
lows that {V^) ^ Tj [aMy^^ <C ("w^) ~ T j {am) and to the leading order we can 
neglect the motions of wormholes which essentially simplifies the above expres- 
sions (e.g., u_ = w = p/ma in ([M]) and dM]))- We also suppose the homoge- 
neous distribution of wormholes in space which gives (rj) = J fw (Xw) dP,^ = 

i?_ — i?+ , 7 ) (here j — byj,U). Then for the scattering term ([50)) we find 



s</™,^-87r / \v\X^{R,p_)-f„,{x,p)]blF{\x~R\,j)dfi, 



where dp, = d^pd^Rd'y. This defines the kinetic coefficient Dm for particles as 
(we use here the property |ti| = |ti'| and d^p — d^p' 

Dm^ — f [W\fL{R,p')-\p\fnAx,P)]9{\x-R\)d'Rd'p, 

ma J 

where g (i?) = / b^^F {R, 7) d'j which gives 

Dm^ — I [K{R + x)-K {x)] g (R) d^i? (44) 
ma J 

with 

Kix,t)= f \p\f{x,p,t)d^p. (45) 



We point out that for the homogeneous distribution of particles K (x) = K 
const and Dm = 0. 

Consider now the second coefficient which is given by 



ma 



(R + x)- (x) 



where 



giR)d^R, (46) 
Ko^ix) = J Po.\p\fix,p)d^p. (47) 
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The multiplier 1/9 appears in ([^^ upon transformations as follows. According 

to dMl) 

J [P» \p'\fix,p')]d'p' = 

= U„p J {p'p~2[p'n)np) \p'\fix,p')d^p'. 

The replacement (|42p means an additional averaging over which gives < 
UaUfj >= ^Sai3, while assuming an isotropic distribution of wormholes over the 
rotation matrix U gives an additional multiplier < Uai^ >= ^Sap- Thus, we see 
that K'^ (x) = {x) in (gS]). The above expressions dH]) and (gB]) show that 
the kinetic coefficients are expressed via the momenta K (x) and (x). 



8.2 functions K{x) and Ka{x) 

Keeping in mind the linearized equations p7l) we will use the following form for 
the quasi-equilibrium distribution function 

n{x,t) I {p- mau{x,t)f \ 

" {2^maT{x,t)f' [ 2rnaT {x,t) ) 

which corresponds to the locally thermodynamic equilibrium distribution and 
n{x,t) = p{x,t)a^ Im. Then we find 

^ = TT— — ^ / exp - "^"^ \p\ d'p. 



(27rmar)3/2 J ^ I 2maT 

Expanding this by u (i.e., K — + K^Ua + ^K'^^UaUp + ...,) we find 

2 , ^ 1 /9 / ma o\ 

K = n— {2maTf'^ [l + —u^) + .... (48) 

To evaluate the function Ka {x) we point out the obvious relations 

d maua -Pa r , . 
^/o {x,p) = /o (x,p) 

which defines the function Ka [x] in the form 
and from we find 



/2mar / 9 maw„\ / ma . 
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Thus, in the leading order by u we get 

K ~ -^(2maT)i/^ Ka{x,t) ~ ^ K {x , t) mau"" {x , t) . 

VTT 3 

Using the relations 

p{x,t) 'mP'a? j fd?p rn^a^ 
01 T — maP I we may rewrite the function K in the equivalent form 



K (x, t) ~ — -j^pcL c. 



2^/2 4 f p \ ^/^ 

PC^s 



For linear perturbations we shall use 5Pm — c^Spm — clpb(t)Sm which gives 
SKix,t)^lK{t)(^c^, + l)^-P, (49) 



and, therefore, we get the functions we are looking for in the form 



^K{t)(l + cl^'^Six,t), SKa,^^K{t)mau^ix,t). (50) 
8.3 Kinetic coefficients Dm and Qm 

In what follows we shall use the parameter Pm/ Pmcl — 1 which gives 5K ~ 
K{t)S{x,t) and K ~ ^^pa^Cs- Then by the use of the above formulas we 
finally find from (|44|) the kinetic coefficients as 



Dm [x, t) = K / [5m {x + R,t)- 6m {x, t)] g (R) d^i? (51) 

ma J 

and analogously from (j46|) 

^ g{R)SR. (52) 



iv„<(a; + i?,0- V„<(a:,t) 



In the last equation the term VqU^Jt, can be expressed from (|34)) as 

^Q^m = -a 3T^m ^i*™ , 

\ at pbQ'^ J 

while for the quantity -j^Dm we find 



m 
pbd^ 



Dm ^ 16V2^Cs J [5{x + R,t)~S {x, t)] g (R) d^R. 
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Consider now some qualitative estimates. Let / g (R) (PR =< &^ > (n„ 
is the density of wormholes in the commoving reference system, the physical 
density is riw/a^). Then we may estimate 

-L)jn ^ 16v27r —d Vyjd 

where ^ (csfe^JT-m) ^, d is the characteristic distance between throats, and 
L is the characteristic scale of the inhomogeneity. For estimates we also can 
suppose Hyj — p^a? /Myj, where — pb^na^b'^ which gives 

Mw _ 4 huj_Pb_ 

Tw ^ 7 2 3 — 

^s^wPwO' o Cs Pw 

and finally we define the estimate for the collision frequency 

TTy/TT bw Pb 

It can be seen that for kcg — > this correction is negligible (the collision fre- 
quency — >■ 0), while for rather short wave-length this term may give the 
leading contribution. 



8.4 Kinetic coefRcients and 

As it was already pointed out the collisions conserve the local number of throats 
which immediately gives the value 

= 0. (53) 

Some non- vanishing value of I?^, < may appear when we take into account for 
the processes of throats merging (annihilation of wormholes) . Such processes are 
important on the very early stages of the evolution of the Universe, since they 
are responsible for the formation of the background distribution of wormholes in 
space. E.g., the origin of the Tully-Fisher relation for spirals was suggested by 
[13] which requires the decay for some portion of primordial wormholes. Such 
a decay is accompanied with an essential reheating, since every merging of a 
wormhole throats radiates the energy ~ M^c^. We however leave this problem 
aside for the future research. 

In the case of the isotropic background the second coefficient does not 
require a separate evaluation. In spite of the fact that during the scattering the 
momentum does not conserves (|20p . some kind of the conservation law for mean 
values takes however place. Indeed, due to the obvious symmetries S\ 4-> S'i 
and t ^ —t 

w (r^, r^, r^, r^) = w (r^, r^, r^, r^) 

we find that ^ 

Qm,a = iPa) = 7^ {Pa ^ Pa) ) 
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where {pa) = J PaStfm,w'Pp- And analogous expression holds for throats 
Qw,a = {P±) = \ {P±-P±) = / P±aStfu],md^P±. Due to the symmetry 
S^l — we can always change (P+) = {P-)- Then using the relation (^0]) 
we find 

l{p-p') + l{P--Pi) = l{il-U)p) + l{il-U)P^). 

For the isotropic distribution the averaging over the matrix U gives < 11^^ >= 
|(5a^ and, therefore, the above expression transforms to — —Qw + + 



iQwt i-e. 



Qni — Q'b 



(54) 



The last equation means not more than the conservation of the total mean 
momentum density. It expresses the balance for the transformation of the mo- 
mentum between particles and wormholes. 



9 The behavior of Unear perturbations 

Consider now the Fourier transform for the density perturbations 
6 A (x) = (27r)- 

Then we find 



-3/2 / ^ikSg^ _^^3,^^ 



:Da {k) = -UA{k)5.r 



where 

^„(fc) = (87r)'c, g{Q)-g[-k] , v^{k)^Q 
Analogously we find 



where 



4 2 



.9(0)-^5(-fc) 



(55) 
(56) 



^Qm (fc) = -^Q- (fc) = ^'^ (I + ^™ (^)) (57) 



(58) 



Since we assume that M^, ^ m, with a very good approximation we may set 
5Pw = 0. Therefore, the master equations are 

I + 2ff + (I + (fc)) <5„ + ^^S^j: + ^^^c = (59) 



I + 1^.,, - (I + -™ (fc)) ^n,, + 
and the Poisson equation (for the Newton potential) 
1 



= 



e-^u = 4vrG {1 + B (fc)) ( + ) , 



(60) 



(61) 



20 



where the function B (fc) was determined by [TT^ and has the forn|f| 



B{k) 



8^(gi(fc)-gi (0)) 
fc2 



(62) 



where gi{X) = f byjF{R,j)dj (so that J gi (x) (fix =< hyj > n^). We point out 
that all the coefficients 17^, z^,„ (fc), and B (k) are functions of the background 
(quasi-stationary) distribution of wormholes F (|i?+ — R- \ , 7) whose exact form 
requires an independent further investigation. 

9.1 Perturbations in the density of wormholes 

Consider first the case when = 0. Then the equation with (pT|) taken 
into account reads 



where Gk = G{l + B{k)). We point out that this equation coincides with 
the standard equation for perturbations in cold dark matter particles. Some 
difference appears however due to the presence of the bias B{k) which reflects 
the polarizability of space filled with the gas of wormholes. Formally such a 
polarizability looks like a scale-dependent renormalization of the gravitational 
constant. By other words, at very large scales (and on the linear stage of the 
development of perturbations) the gas of wormholes reproduces exactly the dark 
matter particles. 

The essential difference appears however on the non-linear stage due to the 
existence of the mutual exchange with the momentum between wormholes and 
baryons (|54p (|57p . It is clear that such an exchange cure the basic failure of 
CDM (i.e., the presence of cusps in the center of galaxies), since dark halos 
around galaxies should rotate. 

9.2 Matter perturbations 

If we set 6 r = 0, then the master equation for perturbations becomes 



The additional coefficients Vm (k) and Qk have the clear physical interpretation. 
Vm (fc) is the collision frequency which describes the processes of the absorption 
and re-radiation of particles by wormholes. It is clear that such processes some- 
what smooth the initial inhomogeneity in the particle number density which 
results in the specific damping we already discussed in [15] . At very large dis- 
tances k ^ d (where d is the characteristic distance between throats) this kind 

^We recall that the bias B consists of the two terms ^ and mUl l. For perturbations the 
leading contribution comes from Bi , while the first part Bq is accounted for by the insertion 
of the rest masses of wormholes. 





AnCkPh S = 0. 



(63) 
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of damping vanishes I'm (k) 0. The additional friction term, proportional to 
flk, describes the mutual interchange with momentum between particles and 
wormholes. We point out that in the linear approximation there exists only a 
transport of the momentum density from particles to wormholes which also leads 
to an additional damping of inhomogeneities. Moreover, this kind of damping 
retains in the long- wave limit as well (f^o 7^ 0). 
Consider the redefinition 



(5g = cxp - / 1/ (k) dt 



then, for the new quantity Sj: we find the standard-type equation 



—5, + 2Hk — St 



9f2 



dt 



(64) 



where 



2Hk = 2H + nk-v [k) , and Gfe = G (1 + E (fc)) . 



By other words, this equation formally coincides with the standard one with 
the renormalized friction Hk and the gravitational Gk constants. The presence 
of the additional damping v (k) and the friction flk somewhat weakens the 
standard Jeans instability. Indeed, in the absence of the expansion {H = 0) 
([64]) defines the dispersion relations {6^ ^ exp (Xt)) 



{X + iyf + 2Hk{X + :^) + 



2^2 



k'c 



AnGkPb 



= 



which defines the two (decaying and increasing) modes S^' ~ exp {Xi,2t) as 



Xi,2 = 7. ± \ 47rGfcPb -I- - 



2 y -r. ^ 

In the limit fc — > we find for the rate of the growth 



Ai 



47rGoPb + ^• 



9.3 Estimates for the background distribution of worm- 
holes 

As it follows from ([55)1 - ([57| . ([5^ the functions v (fc), fifc, and B (fc) depend on 
the background distribution of wormholes in space. The background distribution 
requires an independent consideration and we present it elsewhere, while in this 
section we present some qualitative consideration. 

As it can be seen from (|3ip the function obeys the equation similar to 
the two-point correlation function. We consider here the two different limits 
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(the correct answer lies somewhere in the middle). First, we consider the frozen 
wormholes (i.e., V± = 0). Then the function i^^ — F {\R- — R+\ ,bw,U) is 
completely determined by the process of the wormhole production during the 
quantum stage (presumably an inlationary) of the evolution of the Universe. 
Let us choose the simplest case when ^ S {\R- — R+\ — dw) (e.g., see [TT]V 
Then we get 

g{R)^^^5{R-d^), g,iR) = ^£^S{R-d^) 

where n^jis the density of throats, is the mean throat radius, and d^ is the dis- 
tance between throats. Therefore, we find g (fc) = b'^wnw{2n)~^^^ am{kdw) / (kdw) 
which defines the functions v (k), Qk, and B (fc) in the form 

i.(fc)^-2n.542.)-/^^(l-^), 
^^(k) = (8^)^c.P„n„(27r)-3/2 sin(M^) 



= ^(87r)'c,52^n^(2^)-3/2 



kdyj 
1 sin {kdu 



9 kd,, 



which for kd^ <C 1 give 



"-(2#^6^-(^- 20 

lym (k) « I^c,F^n^(2^)-3/2fc2d^, 



o ^4(8^f 



Thus, we see that in the long- wave limit kd^ <C 1 we have the behavior ^ k'^, 
while B const. 

We point out that this case can presumably be far from reality, since the 
presence of hot matter makes wormholes move chaotically. What we should ex- 
pect that Fyj {V) ~ exp (— M^F2/2r) has the Maxwell-like form. However, due 
to the enormous typical value of the rest mass ^ lO^^rrip, the approximation 
of frozen wormholes can in turn work sufficiently well. 

The second case corresponds to the opposite situation, when all correlations 
between conjugated throats are lost, i.e., when Fw = F {\R- — R+\) ~ const. 
Then we find g (k) ~ go6 (k) and all quantities ly (k), flk, and B (fc) tend to their 
asymptotic (as — ?> 0) values. 
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10 Summary 



In the present section we collect basic results. First of all we showed that the 
background density generates a non- vanishing rest mass of a wormhole (|14[) . 
Therefore, wormholes may play the role of dark matter particles in CDM mod- 
els. Then, we considered the scattering between a particle and a wormhole 
which straightforwardly shows that wormholes move in space. We derived basic 
equations ([T5)) - p^ in the non-relativistic case {V ^ c) which however admit 
straightforward generalization to the relativistic case. Based on the scatter- 
ing equations we suggested the kinetic description in Sec. |6] and derived the 
basic fluid equations in Sec. [71 All those equations are generalized straight- 
forwardly to the relativistic case as well and we present it elsewhere. In Sec. 
[S] we derived equations for linear density perturbations and have shown that 
at large scales density perturbations, related to wormholes, behave exactly like 
standard dark matter particles. However, even on the linear stage there always 
exists the transport of the momentum density from baryons to wormholes. Al- 
ready such a phenomenon is enough to cure the basic failure of CDM particles, 
namely, to remove the cusps in centers of galaxies. Thus, we suppose that at 
present wormholes represent the best candidate for dark matter particles. We 
also demonstrated that development of perturbations in baryons possesses an 
additional damping (with respect to the standard Newtonian instability) . How- 
ever the complete analysis requires considering the relativistic case which we 
present elsewhere. 

As far as basic equations are concern, the generalization to the relativistic 
(GR) case of our basic results seem to have not essential difficulties (Boltzmann 
, fluid equations). The basic difficulty appears when considering the general- 
ization for the bias operator B {x,x') which reflects the polarizability of space. 
For scalar wave equation in the geometric optics approximation such bias was 
considered recently by [16) . However such an approximation works only on 
sufficiently late stage of the evolution of the Universe (where the Newtonian 
consideration works well), while for very early stage it is not sufficient. This 
poses a rather serious problem for future research. 

We point also out that for small astrophysical objects (e.g., jets, nuclei of 
galaxies, etc.) one probably has to account for the dynamics of the additional 
parameters of wormholes (6^, U) which should essentially complexify the theory 
presented. 

The gas of wormholes represents an extremely reach (by physical effects) 
but complex medium. Due to the existence of the polarizability (which has the 
topological origin) it has all properties of a gravitational plasma. It possesses 
also numerous additional non-trivial phenomena and we think it worth saying 
that we have deal with a new state of matter. 
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